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GRAPHICS  AND  COMPOSITE  MATERIAL  COMPUTER 
PROGRAM  FOR  USE  WITH  SPAR 

BY  Gary  L.  Farley  and  Donald  J.  Baker 
Structures  Laboratory,  USARTL 

SUMMARY 

The  SPAR  computer  software  system  is  used  for  finite 
element  structural  and  thermal  analysis.  This  report 
contains  user  documentation  of  additional  computer  programs 
that  have  been  developed  for  use  in  conjunction  with  SPAR. 
These  programs  plot  digital  data,  simplify  input  for 
composite  material  section  properties  and  compute  lamina 
stresses  and  strains.  Sample  problems  are  presented 
including  execution  procedures,  program  input  and  tabulated 
and  graphical  output. 


1 


INTRODUCTION 


Studies  of  structural  configurations  using  finite  element 
models  yield  large  amounts  of  data  which  must  be  analyzed. 
Effective  evaluation  of  these  data  can  be  enhanced  by  a 
graphical  representation.  This  paper  contains  the  user 
documentation  for  computer  programs  developed  or  modified 
to  interface  with  the  data  base  of  the  SPAR  level  14  finite 
element  computer  code  (reference  1) .  These  capabilities  are 
used  in  the  interpretation  of  results  and  in  reducing  input 
to  the  composite  material  section  properties  of  the  SPAR 
computer  code.  These  capabilities  include;  (1)  a  hidden 
line  graphics  program  for  plotting  the  deformed  and  undeformed 
finite  element  model;  (2)  a  contour  plotting  program;  and  (3)  a 
capability  for  the  development  of  composite  material  section 
properties  and  subsequent  postprocessing  of  lamina  stresses 
and  strains  into  a  more  convenient  form.  These  programs 
are  written  in  FORTRAN  IV  language  for  the  Control  Data 
Cyber  series  digital  computers  with  the  Network  Operating 
System  (NOS) .  The  plotting  program  contains  adequate 
comment  statements  to  allow  conversion  to  any  plotting  system. 
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SYMBOLS 


x  Y  2  “  Coordinate  system  fixed  in  model 

*0  Y0  ^0  ~  Coordinate  system  containing  viewing  planes 

z  -  Coordinate  through  the  thickness 


o  a  t 
x  y  xy 


Normal  and  shear  stresses  in  the  elemental 
reference  frame 


°1  °2  T12  -  Normal  and  shear  stresses  in  the  principal 

material  coordinate  systems 

ex  £y  Yxy  ~  N°rmal  and  shear  strains  in  the  elemental 

reference  frame 


Y 
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Normal  and  shear  strains  in  the  principal 
material  coordinate  system 


Ell  E22  E33  -  Three  dimensional  extensional 
G12  G12  G23  and  shear  moduli  and  Poisson's 
UN  12  UNI  3  UN 2 3  ratio 


RHO  -  Material  density 

t  -  Finite  element  thickness 

-  Distance  from  neutral  axis 
k  -  Layer  numbers 
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Program  Capabilities 

The  three  capabilities  described  in  this  report  were 
developed  as  preprocessors  or  post  processors  to  SPAR,  interacting 
only  with  the  SPARLA  data  base.  By  implementing  these 
capabilities  in  this  fashion  the  impact  on  subsequent 
versions  of  SPAR  would  be  minimized.  The  necessary  I/O 
and  data  handling  routines  employed  in  these  capabilities 
are  described  in  reference  (2) . 

Another  common  feature  of  these  capabilities  is  dynamic 
addressing  of  storage  (DAS) .  By  utilizing  DAS  all  problem 
dependent  vectors  of  data  are  stored  in  a  single  working 
vector  in  blank  COMMON.  This  allows  the  user  to  specify 
only  the  central  memory  necessary  to  solve  the  problem.  The 
working  vector,  blank  COMMON,  begins  at  the  first  word 
address  following  the  loaded  program  and  extends  to  the  end 
of  the  available  central  memory  as  defined  by  the  user. 

The  concept  of  element  groups  and  types ,  in  SPAR,  can 
enable  a  user  to  greatly  reduce  the  complexity  of  modeling 
and  interpreting  the  results  of  a  structure.  In  a  similar 
fashion  the  three  capabilities  utilize  this  concept  to  allow 
the  user  to  specify  parts  of  the  structure  to  be  operated 
upon.  This  is  performed  by  the  selection  of  sets  of 
element  groups  and  types.  No  option  exists  to  select 
specific  elements  within  a  particular  group. 
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HIDDEN  LINE  REMOVAL  GRAPHICS  PROGRAM  -  HIDLIN 


The  hidden-line  graphics  program,  denoted  as  HIDLIN, 
is  a  modification  of  that  presented  in  references  3  and  4 
which  was  an  option  to  the  general  orthographic  plotting 
program  of  reference  5.  Plots  of  the  deformed  and  undeformed 
finite  element  model  in  a  3— D  rectangular  Cartesian 
coordinate  system  are  generated  on  a  2-D  viewing  plane  by 
HIDLIN.  The  deformed  plots  are  of  nodal  translations  such 
as  static  displacements  and  vibration  or  buckling  modes. 

Such  plots  are  very  useful  in  debugging  complex  finite 
element  models  and  in  visualizing  the  overall  structural 
response  of  the  model. 

Comparison  of  figures  1  and  2  illustrates  the  clarity 
of  a  composite  cross  beam  model  drawn  by  HIDLIN  versus  a 
drawing  of  all  elements  by  SPAR.  The  deformed  plot  is  shown 
in  figure  3. 

Another  example  of  the  usefulness  of  HIDLIN  is  depicted 
in  figures  4  through  7  showing  the  finite  element  model 
used  in  supersonic  cruise  aircraft  research  at  NASA  — 

Langley  Research  Center  (ref.  6) .  This  model  consists  of 
rods,  shear  webs  and  triangular  and  quadrilateral  membrane 
and  aleotropic  elements.  The  full  undeformed  finite  element 
model  of  the  aircraft  model  is  presented  in  figure  4  with 
figures  5  and  6  representing  the  undeformed  and  deformed  HIDLIN 
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plots  respectively.  Figure  7  represents  the  HIDLIN 
drawing  of  just  the  shear  webs  of  the  model  depicting  the 
internal  structure.  Upon  inspection  of  figures  1  through  7 
the  advantages  of  using  HIDLIN  is  evident  in  debugging  and 
visualizing  the  structural  response  of  the  finite  element 
models. 

Plots  can  be  generated  for  structural  models  containing 
any  combination  of  1,  2  and  3-D  elements.  Faces  of  the 
solid  elements  (3-D)  are  internally  converted  by  the 
program  to  triangular  or  quadrilateral  (2-D)  elements  for 
computation  purposes.  In  terms  of  SPAR  nomenclature  the 
elements  the  user  can  specify  are;  E21-E24,  E31-E33, 

E41-E44 ,  S41 ,  S61,  S81. 

Graphic  errors  in  the  form  of  partially  drawn  hidden 
elements"  can  occur  during  a  normal  execution  of  HIDLIN. 
These  errors,  if  they  occur,  do  not  greatly  detract  from 
the  overall  appearance  of  the  structure  and  do  not  reduce 
the  program's  effectiveness  in  debugging  complex  finite 
element  models  or  in  depicting  overall  structural  response 
as  seen  by  comparing  figures  3  and  8.  To  eliminate  these 
errors,  the  user  can  rotate  the  model  a  few  degrees  or 
adjust  DMAG.  DMAG,  used  in  checking  an  element's  visibility, 
is  the  parameter  that  controls  the  amount  an  element  is 
reduced  about  its  center.  The  causes  of  these  errors  are 
numerical  roundoff  and  a  limitation  on  the  number  of 
segments  of  a  partially  hidden  element. 
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The  additional  computation  required  to  reduce  or  eliminate 
these  errors  is  not  justifiable. 

Input  data.-  The  input  data  deck  is  shown  schematically 
in  Figure  9  and  is  described  in  detail  in  this  section. 

SPAR  data  base  SPARLA  must  be  disc  resident  prior  to  the 
execution  of  HIDLIN. 

SPARLA  data.-  The  SPARLA  data  base  must  contain  the 
basic  structure  topology  prior  to  execution  of  HIDLIN. 


SPARLA  DATA 
SET  NAME 


JLOC  BTAB  2  5 


DEFO  POSI  MASK  MASK 


DEF  E21  MASK  MASK 

"  E24  "  " 

"  E31  "  " 

"  E33  "  " 

"  E41 


"  E44 

"  S41 

"  S61 

"  S81 


DESCRIPTION 

Data  set  containing  nodal 
coordinates 

Data  set  containing  nodal 
translations  in  similar  format 
as  JLOC  BTAB.  This  data  set 
is  used  in  deformed  plots 
only.  This  data  set  can  be 
developed  through  use  of  the 
TRAN  function  in  processor 
AUS. 

Data  sets  containing  element 
connectivities . 
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User  defined  data.-  User  defined  data  includes  a  NAMELIST 
statement  and  data  defining  element  type  and  group  numbers  of  the 
plotted  structure. 

NAMELIST  MAX.-  This  NAMELIST  contains  values  to  allocate 


storage  and  values  specifying  various  program  options. 


FORTRAN  DEFAULT 

NAME  VALUE 


IDISP 


0 


NCKELE 


DSCALE  1.0 


DELX,  DELY  1.0 


KHORZ 


1 


DESCRIPTION 

Deformed  plot  parameter 
0  -  unde formed  plot 
1  -  deformed  plot 
Estimate  of  the  total  number 
of  triangular  and  quadri¬ 
lateral  finite  elements  (must 
be  equal  to  or  greater  than 
the  actual  numbers) 
Displacement  magnification 
factor  used  when  IDISP=1 
Origin  shift  factor  in  the 
scaled  and  rotated  coordinate 
systems 

Integer  designating  the 
horizontal  axis  of  the 
viewing  plane  where  1 =XQ; 


2=Yo'  3=zo 

KVERT  2  Integer  designating  the 

vertical  axis  of  the  viewing 
plane  where  l=Xo;  2=Yq;  3=Zq 
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FORTRAN 

NAME 


DEFAULT 

VALUE 


DESCRIPTION 


PHI  0.0  Angular  rotation  of  model  about 

its  X-axis  in  degrees  (per¬ 
formed  third) 

THETA  0.0  Angular  rotation  of  model 

about  Y-axis  in  degrees 
(performed  second) 

PSI  0.0  Angular  rotation  of  model 

about  its  Z-axis  in  degrees 

(performed  first) 

PSCALE  1.0  Joint  coordinate  magnification 

factor 

DMAG  .99  Reduction  factor  used  in 

reducing  the  size  of  each 
element  about  its  center  for 
checking  visibility  of  an 
element. 

The  scaling  of  the  joint  coordinates  (XYZ)  and  deformations 
(DISP)  and  the  translation  of  the  plotting  origin  (DELX, 

DELY)  is  described  by  the  following  equation; 
original 

XYZpl0t  =  XYZ  model  *  PSCALE  +  DISP  *  DSC ALE  +  DEL (X  or  Y) . 

The  following  card(s)  determine  element  type  and  group 

number  to  be  considered  by  the  program. 
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COLUMN 


FORTRAN  VARIABLE 


1-5  (Right  adjusted)  NGRP  Element  group  number 


NELT 

Spar  Element  Type 

21 

E21 

22 

E22 

• 

L 

E44 

441 

S41 

661 

S61 

881 

S  81 

An  estimate  of  the  field  length  required  to  run  HIDLIN 

is  given  as  follows : 

PL  =  18355  +  6*  (NNOD  +  NCKELE) 
where  NNOD  is  the  number  of  nodes  and  NCKELE  is  the 
sum  of  the  triangular  and  quadrilateral  finite  elements. 

output.-  The  input  NAMELIST  MAX  is  printed  to  verify 
input  data,  followed  by  the  length  of  blank  common  required 
for  program  execution.  The  third  section  of  output  is  a 
listing  of  element  types  and  group  numbers  that  are  being 
plotted.  An  example  of  input  and  output  (printed  and 
plotted)  is  presented  in  sample  problem  1. 

CONTOUR  PLOTTING 


The  contour  plotting  capability  consists  of  two  programs 
which  are  executed  sequentially.  The  first  program, 
denoted  as  STR,  extracts  user  designated  topology  and  stress 
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information  from  the  SPARLA  library  and  sets  up  two  input 

files  for  the  contour  plotting  program,  STCR.  Using  the 

input  data  developed  by  STR,  STCR  plots  contours  over  the 

specified  structure.  Besides  the  contour  levels  only  the 

border  of  the  specified  structure  is  drawn  to  reduce 

confusion.  Three  plots  are  generated  corresponding  to 

a  ,  o  ,  and  t  stresses.  Figure  10  depicts  stress 
x  y  xy  3  r 

contours  of  the  upper  flange  of  the  composite  cross  beam 
structure  shown  in  figures  1  through  3. 

PREPROCESSOR  FOR  CONTOUR  PLOTTING  -  STR 

STR  extracts  stress  and  topology  data  resident  in  the 
SPARLA  library  and  sets  up  two  input  files  for  the  contour 
plotting  program  STCR.  The  user  designates  the  element  types 
and  group  numbers  in  specifying  the  desired  structure. 
Triangular  and  quadrilateral  membrane  and  aleotropic  elements 
are  considered  in  SPAR  nomenclature  as  E31,  E33,  E41,  E43 
elements.  All  the  elements  in  the  specific  groups  are 
included  in  the  input  file  for  STCR. 

Stresses  as  calculated  in  SPAR  are  oriented  in  the  local 
elemental  reference  frame.  These  local  reference  frames 
can  vary  from  element  to  element  and  must  be  transformed  to 
a  common  reference  frame.  This  transformation  is  accomplished 
by  rotating  these  local  stresses  by  the  angle  6  which  the 
user  must  input  into  the  SPAR  material  property  (MATC) 
table . 


11 


Input  data.-  The  input  data  deck  is  shown  schematically 
in  Figure  9  and  is  described  in  detail  in  this  section.  SPAR 
data  base  SPARLA  must  be  disk  resident  prior  to  the  execution 
of  STR. 

SPARLA  data.-  The  SPARLA  data  base  must  contain  the 
basic  structure  topology  and  stress  information  prior  to 
execution  of  STR. 


SPARLA  DATA 
SET  NAME 

JLOC  BTAB  2  5 

MATC  BTAB  MASK  MASK 

SA  BTAB  MASK  MASK 


STRS  E31  MASK  MASK 
"  E33  "  " 

"  E41  "  " 

"  E43  "  " 


DESCRIPTION 

Data  set  containing  nodal 
coordinates 

Data  set  containing  material 
property  table 
Data  set  containing  section 
property  table 

Data  sets  containing  element 
connectivities  and  stress 


information 

User  defined  data.-  User  defined  data  includes  a  NAMELIST 
stateinent  and  data  defining  element  type  and  group  numbers  of 
the  plotted  structure. 

NAMELIST  MAN.-  This  NAMELIST  contains  topology  and  stress 
parameters  pertinent  to  the  execution  of  STR.  MAN  also 
contains  the  appropriate  NAMELIST  parameters  used  in 


NAMELIST  MAX  of  program  STCR. 
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FORTRAN 

NAME 


III 


JJJ 


IPOS 


ICEN 


DEFAULT 

VALUE 


1 


2 


0 


0 


DESCRIPTION 

Horizontal  axis  on  viewing  plane 
where  1=X ,  2=Y ,  3=Z  are  the 
model  coordinate  system 
Vertical  axis  on  viewing  plane 
where  1=X,  2=Y,  3=Z  are  the 
model  coordinate  system 
Location  of  stress  calculation 
for  SPAR  finite  elements 
0  z  =  0  (mid-plane) 

6  z  =  t/2  (upper- surface) 

12  z  =  -t/2  (lower-surface) 

Location  of  stress  component 
(this  parameter  common  to  STR 


and  STCR) 

0  centroidal  stress 

1  nodal  stress 

The  following  parameters  are  required  in  this  NAMELIST 


and  are  passed  to  program  STCR. 


SCLX,  SCLY 
XSHFT ,  YSHFT 
NCONT 
I  LAB 


1.0 

1.0 

5 

0 


HGHT 


.1 


Joint  coordinate  scale  factor 
Origin  shift  factor 
Number  of  contour  levels 
Contour  labeling  parameter 
0  -  No 
1  -  Yes 

Size  of  contour  label 
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FORTRAN 

NAME 


I  COPT 


DEFAULT 

VALUE 


0 


DESCRIPTION 

Contour  specification  parameter 
0  -  program  specifies  contour 
based  upon  the  formula 

SS(J)  =  (NCONT+T  *  (SMAX ( J) -SMIN ( J) ) 
+  SMIN ( J) 

1=1,  NCONT ;  J=l,  2,  3 


Where  SS  is  the  stress  contour 

level  and  J=l,  2,  3  represents 

the  stress  component  o  ,  a  ,  x  , 

a  y  ^y 

respectively 

1  -  user  specifies  contour 
levels  based  upon  the  formula 
SS  (J)  =  (1-1)  *  DSIG(J)  +  SIG(J) 


SIG(J) 

DSIG(J) 

RNDOFF 


0.0  Starting  stress  contour  level 

O.o  Stress  contour  level  increment 

1.0E-7  Roundoff  error  parameter 

used  to  eliminate  small 
deviations  from  0.0  level  contour 
0  Leroy  plotting  option 

0  -  No 
1  -  Yes 

The  following  card(s)  determine  element  type  and  group 
number  to  be  considered  in  this  program. 


I  LEROY 
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COLUMN 


FORTRAN  VARIABLE 


1-5  (Right  adjusted) 

NN 

Element  type 

Spar  Element 

NN 

E31 

31 

E33 

33 

E41 

41 

E43 

43 

6-10  (Right  adjusted) 

NG1 

Element  Group 
Number 

The  user  can  stack  several  different  type  and  group 
cards  to  define  the  appropriate  structure. 

An  estimate  of  the  field  length  required  to  run  STR 
is  given  as  follows: 

FL1q  =  15589  +  3*  NNOD  +  9*  NMAT  +  43*  NSECT 
+  NNN  *  NEL 

where  NNN  =  9  if  ICEN  =  0 

16  1 

NNOD  is  the  number  of  nodes  in  the  model,  NMAT 
and  NSECT  are  the  number  of  entries  in  the  MATC  and  SA 
tables  respectively,  NEL  is  the  number  of  finite  elements 
to  be  plotted  and  ICEN  defines  where  the  stresses  are 
computed  on  the  elements. 

Output . -  The  output  from  this  program  is  in  two  forms, 
printed  and  disk  or  tape  resident.  Input  NAMELIST  MAN  is 
printed  to  verify  input  data,  followed  by  the  length  of  blank 
common  required  for  program  execution.  The  third  section  of 
output  is  a  listing  of  the  element  types  and  group  numbers 
that  are  to  be  plotted.  The  disk  resident  output  from  this 
routine  to  be  used  in  the  contour  plotter  is  formatted  in  card 
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images  and  located  on  tape  9  and  tape  10 .  The  data  structure  of 

tape  9  consists  of;  NAMELIST  MAX,  sequential  ordering  of 

joints  and  element  connectivities.  Tape  10  contains  the 

elemental  centroidal  or  nodal  stresses  (a  ,  a  ,  t  )  . 

x'  y  xy ' 

CONTOUR  PLOTTING  PROGRAM  -  STCR 

The  contour  plotting  program,  STCR,  draws  the  border 

of  the  specified  structure  along  with  specified  contour 

levels.  Those  contours  intersecting  the  border  can  be 

optionally  labeled  corresponding  to  the  printed  output. 

The  contour  levels  are  tabulated  in  the  program  output.  The 

method  of  drawing  contours  employed  in  STCR  is  analogous  to 

that  of  reference  5  and  will  not  be  discussed  here.  Three 

contour  plots,  depicting  a  ,  o  ,  and  x  ,  are  drawn  during 

x  y  xy 

a  program  execution  as  demonstrated  in  sample  problem  2. 

The  contours  depicted  in  figures  10a  through  10c  have 
sharp  corners  and  small  extensions  which  could  mislead  the 
interpretation  of  results  by  a  novice  user  of  a  finite 
element  program.  These  corners  are  a  result  of  the  stress 
averaging  at  the  element  nodes  and  calculation  of  the  contour 
line  segment.  These  contours  can  be  made  smoother  by  using 
a  finer  mesh  of  elements. 

Input  data.-  If  program  STCR  is  executed  sequentially 
after  STR,  no  additional  input  is  necessary.  The  following 
data  are  given  to  allow  this  program  to  be  used  alone.  Tape  9 
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(input  file)  contains  NAMELIST  MAX,  sequential  ordering  of 
joints  and  element  connectivities  while  tape  10  contains 
element  centroidal  or  nodal  stresses  to  ,  a  ,  t  ) . 

x.  y  xy 

NAMELIST  MAX.-  This  NAMELIST  contains  all  the 
parameters  necessary  to  control  the  contour  plotting. 


FORTRAN 

NAME 

DEFAULT 

VALUE 

DESCRIPTION 

NNOD 

Number  of  nodes 

NEL 

Number  of  elements 

See  NAMELIST 

MAX  of  program 

STR  for  description  of 

following:  SCLX,  SCLY,  XSHFT,  YSHFT ,  NCONT,  ILAB,  HGHT, 

ICOPT,  SIG ( 3) ,  DSIG ( 3) ,  RNDOFF,  I LEROY . 


♦Sequential  ordering  of  nodal  coordinates 

X  ,  Y  ;  (NNOD  cards)  (2F10.4) 
c  c 

•Element  Connectivities  (NEL  cards)  (415) 
•Element  centroidal  or  nodal  stresses 
if  ICEN=0  centroidal  stresses 


a  ,  a  ,  t  (3F10.4) 
x'  y  xy 


if  ICEN=1  nodal  stresses 
a  (4F10.4) 

X 

a  (4F10.4) 


x  (4F10.4) 
xy 


An  estimate  of  the  field  length  required  to  run  STCR 


is  given  as  follows : 

FL1q  =  18088  +  6*  NNOD  +  4*  NEL 
where  NNOD  and  NEL  are  the  number  of  nodes  and  finite 
elements  respectively. 
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Output.-  NAMELIST  MAX  is  printed  to  verify  input.  The 
length  of  blank  common  is  indicated.  The  contour  number  and 
its  associated  value  is  listed.  A  typical  plotted  output  is 
shown  in  sample  problem  2. 

COMPOSITE  LAMINATE  CAPABILITY 

The  composite  laminate  program  consists  of  two  programs; 
SAT,  which  is  used  to  calculate  the  section  property  table 
(SA  data  set) ,  and  program  STST  for  the  determination  of 
lamina  stresses  and  strains  in  the  elemental  and  principal 
material  directions.  The  input  of  composite  laminate 
section  properties  into  SPAR  is  possible  by  direct  input  of 
the  laminate  stiffness  coefficient  matrix  or  by  ply-by-ply 
specification  of  the  lamina  stiffness  matrix.  Either  method 
requires  the  user  to  execute  additional  programs ,  separate 
from  SPAR,  to  generate  this  input  data.  In  the  execution 
of  program  SAT  the  user  defines  the  necessary  orthotropic 
material  properties  and  the  appropriate  stacking  sequence  for 
each  section  while  the  program  performs  all  necessary 
calculations  for  the  SPAR  SA  table  entries. 

The  stress  recovery  in  SPAR  has  two  forms.  For  all 
elements  with  section  types  other  than  LAMINATE  the  laminate 
stress  resultants  or  average  laminate  stresses  are  computed 
while  the  LAMINATE  section  type  provides  for  the  calculation 
of  stresses  on  a  ply-by-ply  basis  only.  Neither  of  these 
allow  for  the  direct  determination  of  strains  in  the  elemental 
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or  principal  material  reference  direction.  In  addition  to 
the  SPAR  generated  stress  data  the  user  can  have  program 
STST  calculate  stresses  and  strains  in  the  elemental  and 
principal  material  directions  on  a  ply-by-ply  basis.  This 
information  can  be  optionally  directed  to  a  file  in  a  format 
consistent  with  input  for  the  contour  plotting  program  STCR. 

SHELL  SECTION  PROPERTIES  PROGRAM  -  SAT 

The  SAT  program  generates  shell  section  properties,  SA 
table  entries,  applicable  to  composite  materials.  Besides 
the  SA  table  data  set,  developed  by  SAT,  a  data  set 
containing  constituent  section  properties,  denoted  as 
ARMY  COMP  (NLAY)  0,  is  stored  in  the  SPARLA  library.  These 
data  sets  in  addition  to  the  stress  resultants  computed  by 
SPAR  are  used  by  program  STST  in  computing  lamina  stresses 
and  strains  of  the  appropriate  elements.  The  SA  tables 
generated  by  SAT  are  applicable  to  the  following  SA  sections; 
isotropic,  membrane,  plate,  uncoupled  and  coupled. 

The  user  specifies  the  constituent  material  properties 
to  be  used  in  developing  all  the  different  section 
properties.  Each  different  section  is  defined  with  respect 
to  ply  orientation,  lamina  thickness  and  constituent 
materials.  Sign  convention  and  consistency  with  the 
mathematical  formulation  of  SPAR  is  maintained.  Figure  11 
depicts  sign  convention  and  typical  laminate  construction. 
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Input  data.-  The  input  data  deck  is  shown  schematically 
in  figure  12  and  is  described  in  detail  in  this  section. 

SPAR  library  SPARLA  must  be  disk  resident  prior  to  the 
execution  of  SAT. 

Comment  card.-  One  card  required  to  identify  the  data 


de  ck . 

NAMELIST  MAX.-  This  NAMELIST  sets  up  the  required 
number  of  cards  to  be  read  later  for  different  section 
properties.  Note,  there  are  no  default  values  in  this 

namelist. 


FORTRAN  VARIABLE 
NSECT 


NLAY 


DESCRIPTION 

Number  of  different  section 
properties  to  be  developed 
Maximum  number  of  layers  in 


NMAT 


NAMELIST  PROP.-  This 
identification  number  and 


a  section 

Number  of  different  material 
properties  to  be  input 
NAMELIST  contains  the  material 
its  associated  properties. 


FORTRAN  VARIABLE 


I  MAT 


Ell,  E22,  E33 , 
G12 ,  G13 ,  G23, 

UN 12 ,  UN 13 ,  UN 2  3 

RHO 


DESCRIPTION 

Material  property  identification 
(number  not  to  exceed  NMAT) 

Three  dimensional  material 

properties 

Material  density 
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NAMELIST  SECT.-  This  NAMELIST  contains  the  section 
topology  information  to  be  read  in  NSECT  times. 


FORTRAN  VARIABLE 

I  LAY 

ISPARM 


DESCRIPTION 

Number  of  layers  in  section 
SPAR  material  number 


to  jviatc  table 


ISPTYP 


ZSHFT 


entry 

SPAR  section  type 
ISPTYP  SPAR  SECTION  TYPE 
membrane 
plate 
uncoupled 
coupled 
isotropic 


1 

2 

3 

4 

3 


Neutral  axis  shift 
The  following  layer  identification  card(s)  to  be  read 
ILAY  times  for  each  section. 


COLUMN  FORTRAN  VARIABLE 

1-10  THETA (I) 


11-20  T ( I ) 

21-25  (Right  M(I) 
adjusted). 


DESCRIPTION 
lamina  rotation  angle , 
relative  to  the  elemental 
and  material  coordinate  systems, 
refer  to  figure  11(b) 
lamina  thickness 
lamina  material  number  (IMAT) 
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SPAR  stress  recovery  cards 


FORTRAN  VARIABLE  DESCRIPTION 

F(I,  J)  SPAR  stress  recovery  parameter 

1=1,  2,  3;  J=l,  2,  3,  4,  5,  6 
(18  entries)  Format  (8E10.3) 

An  estimate  of  the  field  length  required  to  run  SAT  is 
given  as  follows: 

FL1q  =  18867  +  37*  NMAT  +  (12*  NLAY+140)  *  NSECT  +  NLAY 

where  NMAT  is  the  number  of  different  materials,  NLAY 
is  the  maximum  number  of  layers  in  any  laminate  and  NSECT  is 

the  number  of  sections  (SA  table  entries) . 

Output . -  All  input  NAMELISTS  are  printed  along  with  the 
[ABD]  matrices  (a  6  x  6  laminate  coefficient  stiffness  matrix) 
and  their  inverses  for  each  section  (SA  table  entry) .  An 

example  of  input  and  output  is  presented  in  sample  problem  3. 

STRESS-STRAIN  CALCULATION  PROGRAM  -  STST 

Program  STST  determines  lamina  stresses  and  strains  in 
the  elemental  and  principal  material  directions.  A  contour 
plotting  option  is  available  to  the  user  in  the  creation  of 
a  stress  or  strain  file  applicable  to  program  STCR.  As  was 
the  case  in  programs  HIDLIN  and  STR  the  user  defines  the 
structure  of  interest  by  specifying  group  numbers  and  element 
types.  All  the  elements  within  such  a  specification  would 

have  their  lamina  stresses  and  strains  computed.  Successive 
specifications  can  be  defined  to  establish  the  desired 

structure . 
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Input  data.-  The  input  data  deck  is  shown  schematically 
in  Figure  9  and  is  described  in  detail  in  this  section. 

SPAR  library  SPARLA  must  be  disc  resident  prior  to 
execution  with  data  set  ARMY  COMP  (NLAY)  0,  generated  by 
program  SAT,  resident. 

SPARLA  data.-  The  SPARLA  data  base  must  contain  the 
constituent  material  properties  for  each  section,  as 
generated  by  SAT,  and  stress  information  prior  to  the 
execution  of  STST. 


SPARLA  DATA 
SET  NAME 

ARMY  COMP  (NLAY)  0 


STRS  E31  MASK  MASK 
"  E33 

"  E41  "  " 

"  E43  " 

MATC  BTAB  MASK  MASK 


DESCRIPTION 

Data  set  containing  the  consti¬ 
tuent  material  properties  for 
each  section 

Data  set  containing  element 
stress  information 


Data  set  containing  material 
property  table 

User  defined  data.-  User  defined  data  includes  a  NAMELIST 

statement  and  data  defining  element  type  and  group  numbers 
of  the  plotted  structure. 

NAMELIST  MAX.- 


FORTRAN  VARIABLE 

I  PLAY 

IGRPH 


DEFAULT  DESCRIPTION 

1  Layer  number  for  plotting 

0  Plotting  parameter 

No  =  0 


Yes  =  1  =  a  ,  a  ,  x 

x  y  xy 

2  =  e  ,  e  ,  y 

x  y  '  xy 
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FORTRAN  VARIABLE 


DEFAULT 


DESCRIPTION 


.5  Ratio  defining  where  the 

stresses,  though  the  thickness, 
in  each  layer  are  computed. 

PER  is  relative  to  the  lower  edg 
of  each  lamina.  NOTE:  PER  =  .5 
is  mid-surface 


The  following  card(s)  indicate  the  element  type  group 
number  to  be  considered  by  the  program. 

COLUMN  FORTRAN  VARIABLE  DESCRIPTION 


(Right 

NELE 

Element 

type 

adjusted) 

NELE  SPAR  ELEMENT  TYPE 

31 

E31 

33 

E33 

41 

E41 

43 

E43 

i  (Right 
adj  us ted) 

NGRP 

Element 

group  number 

An  estimate  of  the  field  length  required  to  run  STST 
is  given  as  follows: 

FL1q  =  24249  +  24*  NMAT  +  NSECT*  (11*  NLAY  +  44) 
where  NMAT  is  the  number  of  different  materials,  NSECT 
is  the  number  of  sections  (SA  table  entries)  and  NLAY  is 
the  maximum  number  of  layers  in  any  laminate. 


Output.-  The  input  NAMELIST  MAX  is  printed  out  along 
with  information  pertaining  to  laminate  construction 
generated  by  SAT.  For  each  element  of  the  specified  group 
and  type  the  laminate  strains  and  curvatures  and  lamina 
stresses  and  strains  in  the  elemental  and  principal  material 
directions  are  printed.  Stresses  or  strains  for  designated 
layers  are  printed  on  tape  10  for  use  in  conjunction  with 
program  STCR.  An  example  of  input  and  output  is  presented 
in  the  sample  problem  4. 


Concluding  Remarks 

The  computer  codes  described  in  this  report  have 
been  found  to  reduce  hand  manipulation  of  data  and 
improve  visualization  of  results.  The  HIDLIN  program 
proved  beneficial  in  debugging  complex  finite  element 
models  and  visualizing  the  overall  response  of  the 
model.  The  composite  material  programs  SAT  and  STST 
greatly  reduce  the  manipulation  of  input  data  and 
extend  the  computational  capability  of  SPAR.  The 
contour  plotting  program  used  in  conjunction  with 
SPAR  or  program  STST  significantly  adds  to  the  stress 
plotting  capabilities  currently  available  in  SPAR. 
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SAMPLE  PROBLEMS 


Problem  1 

This  example  illustrates  the  input,  output,  and  typical 
plot  from  program  HIDLIN.  The  Calcomp  plotter  with  Leroy 
pen  was  utilized  to  plot  figure  2.  A  listing  of  input  data 
cards  follows: 

>max  PHT=-4Ff  theta* ♦4*5,  fPSi *45 * NCKELE=450, lDlSP=0f 0SCALE=1 00 . f  ILF;*0Y  =  1  SEND 
1  43 

?  4  3  . .  -  .  '  . 

3  4  3 

4  43 

6  *-3 

7  4  3 

P  43 
9  43 

10  43 

16  43  .  _  . . .  —  •  - .  -  - . 

11  43 

1?  43 

13  43 

14  43 

15  43 

17  ..  43  ' .  . . . - . . -  . - - - - - 

1  a  43 
1  33 

?  33 
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The  output  listing  for  Problem  1  follows. 


a  M  A  X 

T0T5P 

i « 

MCKFL.F 

= 

4  5  0  « 

DSCAI  F 

rr 

.  1 F  +  0  ?  « 

OF  lx 

s 

.  1 E  ♦  0  1  t 

OF  I  V 

= 

.  1 F  ♦  0  1.  » 

KH0P7 

= 

1  t 

KVFPT 

?« 

OHT 

-,45F*0? 

ThFTA 

.451:  *0?. 

PST 

= 

.  4  5  E  +  0  ?  « 

P5CAI.F 

- 

.  1  F  +  o  1  . 

PM  a« 

= 

.Q9F  +00  , 

It.FSOY 

= 

1  ♦ 

5EM0 

STOPAGF  PFOUIRFO  IN  COMMON  (DEC.)  5Q10 

STOPAGF  PFPUTFFO  T N  COMMON  (OCT* ) 0000750G 


GPO'IP  NIJMpFF- 

1 

FLFMFNT  TYPF  = 

43 

GPOHP  Ml  IMF  F P  = 

? 

F l  FMFMT  TYPE:  = 

43 

GPO'IP  pi!MPFF  = 

3 

FlfmFMT  TYPF= 

43 

GPOI  IP  Ml  |MRFk  = 

4 

FLFMFMT  T  Y  P  F  = 

4  3 

ppnuP  M 1 1 M p F F  = 

s 

FLFMFMT-  TYPE  = 

4  3 

'=d  m()MRF'P  = 

6 

Of  i 

l 

FLFMFMT  T YPF  = 

3  3 

GPOIJP  Ml  IMP  F  P  = 

? 

Problem  2 


This  example  illustrates  the  input  and  output  data  from 
programs  STR  and  STCR.  A  Calcomp  plotter  with  Leroy  pen  was 
utilized  to  plot  (figure  10)  the  stress  contours  on  the  upper 
flange  of  a  composite  cross  beam.  A  listing  of  input  data 
cards  follows: 

fMEM  TTi  =  i,jjj  =  3,MCONT  =  ?0»ILAfl=l,ICOPT=t,Sir.(l)=-4S0n,Sir.(?)=-400n,ILFRriY  =  l, 

SIP  (3>*-lft00«OSIG{n=500.OSIG<2)=500»OSIG<3)=2<»0»ICENs0,SCLX*.5»SCLY*.5 
4  3  6 

43  7 

43  P 

43  14  . .  . . .  .  . - .  -  . - . . . . . .  . 

33  ? 
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The  output  listing  for  program  STR  of  Sample  Problem  2 
is  shown  below. 

%  I V,  ft  h 

1  J  T  =  1  , 

•JO  J  =  .?♦ 

SCL>  S  .^.F+nOt 

sci  v  -  ,cf+nn. 


>Hrr  =  .1F  +  PJ, 
vc  hft  =  .  IF  *n  J  . 

TC»  =  f i « 

f  rn‘  T  =  POf 

Tcr^  =  c  • 

il.  ,/st-:  =  1. 

K-t-'T  =  .  K  +  oo, 

[C.rpT  r  i. 

''  1  =  -.4CC+04,  -  .  4  F  +  ri  4  ,  -.IF +  04. 

I '  s  T  <"  -  .  i  +  0  3  .  .Fr.no,  .  R  F  +  0  3  , 

TI.FfOY  =  1, 

u- f  1 1  n o F F  =  .  H  -Of. 

C  F>ifi 

STO'-MSF  T*  C0P'*OV  +  7  OP  ]  ft  X  r.jer  | 

cTr-RAf-f  usf  ,o  |w  co,‘-/“'nr.i  +  7  OR  ?n  x  mfi 


F 1 

F 

{;P»- 

•T 

T  VP  h  r 

4  ^ 

r-  i 

■  ri 

IIP 

r  ! 

i m f  n-  = 

h 

F 1 

r 

V;P!' 

T 

T  V  n  f  = 

4  i 

SF 

'•  O 

}  }D 

M 

i  *«  opt-- 

7 

F  1 

F 

m  r 

T 

T  YPF  = 

4  .  -4 

HP 

!  IP 

*  1 

*  p  p  P  = 

H 

F! 

F 

M  FF- 

T 

T  YPF  = 

4  4 

r;  r. 

n 

1  IP 

M 

-  *  P  F  ^  = 

1  4 

F'l 

F 

U  F  K! 

T 

T  Y  P  F  = 

3  P 

(OFT.)  ?31R 

(orr. )  o  u  o  o  4  6 ] 7 
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The  output  listing  for  program  STCR  of  Sample  Problem  2 

is  shown  below. 

"  A  y 

on  =  4  03, 

1  =  70, 

TCF*  =  n, 

>  >  T  r  .  ]  F  +01* 
vnorr  =  .It *01, 

‘‘■Cl*  -  .  4  !■:  +00. 

‘civ  -  .c  f  +  n  o , 

-co  nt  =  ?n. 

1 1  a  0  =  1  , 

h FT  =  .  ]  r  +  0  f !  , 

‘TO  -  -  ,  4 't  !■  +  04,  —  .  4  F  +  ^  4  .  - .  i  p  +  o  4  , 

■  c-  T  3  =  .  r >  +  0  3  ,  .  o  F  +  03.  .  ?  F  +  o  3  , 

T  C  (  r  T  =  1. 

T!  n-ov  =  ]. 

■  ’  r  nor  = 


"  FOi; 

STOPAOf  licpr,  p-  OfiMyriM  (OFC.)  ?7?3 

STOP  a  OF  itSFl  IF;  COM‘«n\'  f  OCT  .  )  0  0  0  0SP47 


1 A  *  -  •  T  >  STPfSF  VAUJFS 

.0470+03  -.7440+04 
.  OP?c  +04  -.F43F+03 

.4740+03  -,F7hf+03 


*=■  s  M 

U  L 

i 

rnMoup 

KJ>-  I 

-00  1 

V  A  1.3  iO 

-.4000000F  +04 
'■  °  o  r  0  0  0  +  0  4 

r  >  ■* 

-  ■ 

'  '  r^r*  S 

V  A  U  1 0 

*  *  * 

f  Sc 

CO  OTO\Jt{ 

Mijfo 

O  E  F?  1 

V  A|J  IF 

-.00000  0  t1  [. 

F  So. 

0 

c.  n^TOiiP 

MUM 

m  pp  4 

VAIJiO 

-.40000000+03 
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Problem  3 


This  example  illustrates  the  input  and  output  from* 
program  SAT  which  computes  the  SA  Table  entries  for  SPARLA. 
A  listing  of  input  data  cards  follows: 


F  CASF  FOP  PROGRAM  SAT 
«.MAY  fcFrT=?»NLfiV  =  ft,MMAT=?  TFNP 

■jppop  IMAT=ltFll=17.E6.E??»2.E6.G12».52£6.UMl*».38  SEND 
«PROd  T  wa  t  =  ?  » F 1  1  =  9 .  EG  «E??=9  .E6  *  G 1 2=  •  6Ef  *  UNI  ?=  •  3  SEND 
<SFCT  ILAY=4f X^PAPMr] , !SPTYP=4  $END 


0.0 

1.0 


45. 

.0455  . 

?  . 

.  ,  - - 

-45. 

#  0 455 

? 

-45. 

.0455 

? 

4*. 

.  0455 

2 

1.0 

1.0 

1.0 

o 

. 

o 

1.0 

0 . 0 

0.0 

0.0 

0.0 

0 . 0 

’  TLAYs5* 

ISPARMs^, 

ISPTYP=4 

?FND 

0.0 

.033 

1 

45. 

.  0  ?6 

2 

o.o 

.044 

1 

A  c  . 

•  0  ?5 

? 

o.O 

.033 

.  1 

_  .  .. 

1  .0 

1  .0 

1  .0 

0.0 

1.0 

0.0 

o 

• 

o 

o.n 

0.0 

0.0 

0.0 

1.0 


0.0 

1.0 


o.n 

1.0 


1.0 

1.0 


1.0 

1.0 


1.0 

0.0 


1.0 

0.0 
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The  output  listing  for  program  SAT  of  Sample  Problem  3 
is  shown  below. 

*  v  \  v 

^frT  =  ?. 

M.f  v  a  <*♦  . . . .  '  . .  "  •  -  ■ 

*•*>  f  7  =5,  .  "  . ” . .  \7"' 

**f  I’&f.P  CO*vfifci  II>EC.)  ...  SOS . 

5T0*>4f.f  I'SFft  TK  COV^:  (OCT  •)  0000ft  771  ..  .  .  ...  7 


* :  r 
.  1 


F!1 

.17  3**0?.  ... 

.  o  ft  ft  F  *  o  7 

l  M  ■> 

*73?’  o  (1  ♦  ft  ,i 

,  ~  ft  o  r  r  >  f  ♦  n  o 


r?? 

F.33 

si? 

•  2CrE*07.  G 

_ *5.^.0  £*ft  S 

.  <5  o  n  fc' 1 0  7  fl  , 

,4(;"r*  ns 

HMP 

i 

ft. 

ft. 

0. 

(i  13 


G?3 


PHO 

0* 

0. 


TLAV  * 

!C^A*V  =  1. 

tSt'TVp  s  4 

n-FT  -  0.0, 


THfT:  ThTC^VPs’S 

-  „  4  *  ft  r  +  c  ,  4  «.  £  c  -  ft  * 
#4SAp*pj?_  ,r;-r.r«f'i 


^AT 

? 


♦  r  r  v  t  r  f 

K  *f  ♦  <! 

r  V 1  W  v  P  ?.  v.‘  A  >  •  t  ^  t  r  y  s  ( F'  j 

.  10 or ♦<'!  *ino->oi  o. 

ft,  0 

.  .loot-:* 

01  .mcF*01 

i  *■  *  *  *  \  n 

«  ft  • 

.  .  1  ft  o  r  *  M 

ni  o. 

if 

.  i  r  *f>? 

i 

. jrAOho  nf  +  07 

ft. 

•  4<i«71F-n»-» 

.4.01*4  01  Of 

MfuSr-?4 

.  10  s'  rtK^ri.  +  ft  7 

.  1?7«;'IK‘F*0T 

0. 

.*gw,Q]  pf  -op 

,r  1 4<yft7.i^-f,g 

*p?71P0SF-?4 

T  * 

#  A  3  i  v  0  *  O  tP  *  0  0 

,H?7J H04P-P V  - 

,P?71 00SK-?4 

-,4ft7AS3sF‘-09 

.*'P»!4'i }  Pf -09 

.0^71 r -?4 

,3S7J  01 «F*04 

*?•»?*- lh?P-*  04 

-,#4?)Rl?K-l 1 

'  P?  «*fto 

,  A  1  4  W n  ,f 

-.P -'71  *fttt  -?A 

.  PUpCll^pr  *r>4 

,  is 7*  01  .(,4 

,8^?! Pi ?f-J 1 

*  "* 

-.4«!74S-»M:  -no 

?r  -i  1 

,M4?lHt?f-l 1 

.  1739ft lftf>  04 

rocp  or 

.  1  ' 

.11'*’ 

#  A}  *.|'J  7P»  -  ^ 

A  V  'O  WAT-  tv  V 
,  v1?.  •■•  7  ■?  ft  w  ,.n< 

•  •  ‘f  T  f  O’1: 

t.>  ^rcT ! 

1 

‘  '  '  r  .O^F.iV 

? 

4  :V»oi  4«>  -?r. 

-  .*Ot  1  S7?fT-l  A 
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Problem  4 


This  example  illustrates  the  input  and  output  of  program 
STST  lamina  stress  and  strains.  This  example  uses  the  file 
ARMY  COMP  ( NLAY )  0  from  problem  3.  A  listing  of  input  data 
cards  follows: 

*  i ..  a  v  v  •  f  c*  =  n  #  n  ft  p  t  ■ ) 

?  ^  1 
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The  output  listing  for 
is  shown  below. 


$MA  X 

I  play 

*  1. 

IG9PH 

*  0* 

PFO 

a  0.0 

5-ENR 

NUMBER  OF  MATFPIALS  a  3 

maxivjm  number  of  layers/ section  *  20 

Nt»vepp  oF  S*  TAHLE  ENTREES  a  23 


S TORAP-E  USED  IN  COMMON  (DEC.)  6145 

STORAGE  USED  IN  COMMON  (OCT  * )  0981*001 


G»OU°  nuvrER  1  ELEMENT  type 


10 

LN 

SOX 

EPX 

SGI  fpi 

SGY 

MIP-PLANE 

STRAINS  .668E-03  -.870f  «>-  3  , 

49AE-93 

X 

l 

-57. 

.00021 

-10777.  -.*0063 

-10777, 

1 

p 

3726. 

.00024 

37?*.  .  O.-0  24 

“  *  1 53  , 

1 

3 

op. 

.00031 

-1159S.  -.A.jr.6 A 

-11 59ft. 

1 

4 

534). 

.00034 

5  3*4*  .  ,»:>034 

“1*6 5. 

1 

«; 

?<M. 

,00017 

-12  1*5.  c  7 2 

-12146. 

1 

* 

6C55. 

,00044 

5955,  4* 

-1  I 96  , 

1 

7 

?r'6. 

.00047 

-1  2  9*  r*  •  -.i‘>i77 

-12966. 

1 

8 

8032. 

.00060 

8  0  3  2  .  •  s*  "i  0  5  0 

-1217. 

1 

Q 

511. 

♦ono5f 

-13785.  -.0^0*2 

-13^86. 

1 

10 

-709. 

.000*0 

212.  , ft. 139 

“2:43!. 

1 

11 

-197, 

.00073 

-2835,  -.n.'iOSA 

-?!63. 

1 

1? 

315. 

.oops* 

-2.'SP.C  —  ,  5 i 

“1896, 

1 

13 

e?6. 

.  0  0  0  C  9 

1CM.  ,r.-039 

“1627. 

MO-PLANE 

STRAINS  .6602 

-03  - . R 7 9 F - *; 3 

4 98 F- 03 

p 

1 

“r,7  • 

.00021 

-1  07*^7  . 

-10777. 

p 

? 

3736. 

,00074 

3735.  .00024 

“1133. 

;> 

3 

ce . 

.0003] 

-H5150.  -.0J3-VI 

-nsof . 

? 

4 

6341. 

.00034 

5341.  .00034 

“  1 1 6  5  • 

2 

6 

201  . 

.00037 

-121—5.  -.00 072 

-12146. 

2 

6 

6°c6  . 

.00044 

69  65.  .76044. 

“1196. 

2 

7 

.  366. 

.00047 

-12965.  -.00077 

-12965, 

A 

003?, 

.  0*0  50 

8032.  .nonso 

“1217. 

611  . 

.00056 

-13785.  -.0008? 

-13785. 

.  - 70c . 

.00060 

-3351,  -.'•.'064 

“2431. 

1* 

’97, 

.03073 

495,  .non 39 

“2163. 

14 

3 

,000*6 

778.  .00039 

-1806, 

1  4 

4 

-1662,  -.00045 

—  1 6?  7  - 

u 

6 

-  ■ 

14 

6 

57**c, 

14 

7 

-7463 . 

—  .00,..  , 

•  4 

0 

61e2. 

-.00022 

1  /  .  ,.*•!!  f  1 

1  ic>?6. 

14 

Q 

-76Qa , 

-.00027 

-11864,  -,00144 

“2870. 

u 

1  0 

134, 

-.00029 

421,  ,00023 

2?9ft. 

1  4 

11 

-irr.o. 

-,00039 

149?,  -.00-0  01 

1709. 

14 

12 

-2663. 

-.00049 

97?.  ,00002 

llio. 

1  4 

13 

-3896, 

-.00059 

-3819.  -.00031 

530. 

M!0-PLANF 

STRATUS  - , 3  4  0 1 

-03  ,2646-03  . 

-n? 

16 

1  -169S6. 

-.00012 

-2*988.  -*00333 

“Tsoi . 

16 

? 

9  ?  5  3  , 

-  .  0  0  0  l  4 

30860,  .0835] 

1974?, 

16 

3 

16260. 

-.00017 

-26869.  -.00315 

“7  5  4  0. 

15 

4 

544?  , 

-.00018 

28471.  #so 32* 

18 199, 

15 

A 

I 664?. 

-  .  0  0  n  ?  p 

-25123,  -.10303 

“7067, 

16 

6 

763?. 

-.00023 

26082.  .0030! 

1 6*6*  , 

15 

7 

16215. 

-.00025 

-24004.  -.Or? PS 

“6,806. 

15 

ft 

7091  . 

—  ,  0  0  0  ?  4 

24489.  ,00284 

1 6*?p. 

15 

9 

1 4  5  c’  0  . 

-.00029 

-22885.  -.00-67 

“*546. 

1  a 

10 

460  , 

-.00031 

-1381.  -.rro 64 

“76?. 

l  5 

1 1 

-547. 

-.'0  Oft  37 

-7.  .00045 

“  c  9  2  . 

I  6 

1? 

-1551. 

- .  0  ft  0  4  3 

-110?.  .00028 

- 1  ?  3  ?  . 

16 

13 

-?659. 

-.00049 

-1835,  —  •  0  0  4  2 

“ 1 47  ?  . 

program  STST  of  Sample  Problem  4 


EPY 

SG? 

EP2 

TAUXY 

GMXY 

TAU12 

CURVATURES 

.?90E-0? 

1 157E-02  - 

.678E-03 

-.00063 

-57. 

.00021 

314. 

.00060 

-314. 

-.00066 

-1133. 

-.00065 

310. 

.00060 

310. 

-.00068 

98. 

.00031 

303. 

.00058 

-303. 

-.00070 

-1166, 

-.00070 

299. 

.00057 

299. 

-.0007? 

20). 

.00037 

295. 

,00057 

-?95, 

-.00076 

-1196. 

-.00075 

287. 

.00055 

287. 

-.00077 

366. 

.00047 

283. 

.00054 

-283. 

-.00079 

-1217. 

-.00079 

279. 

.00054 

279. 

-.0008? 

5)1. 

.00056 

27?. 

.0005? 

-272. 

-.00084 

-3351, 

-.00064 

1782. 

.00051 

-861  . 

-.00091 

495. 

.00039 

1675. 

.00048 

983. 

-.00098 

778. 

.00039 

1568. 

.00045 

1105. 

-.00106 

-1862. 

-.00045 

1461. 

.0004? 

-1227, 

CURVATURES 

,?9or-r? 

157E-0? 

.67BF-03 

-  .000*3 

-5  7. 

.00021 

-314. 

-.00060 

314. 

-.000*5 

-1133. 

-.00066 

-310. 

-.00060 

-310, 

-.000*8 

98. 

.00031 

-30.3. 

-.00058 

303. 

-.00070 

-1165. 

-.00070 

-299. 

-.00057 

-299. 

-.00072 

201. 

.00037 

-295. 

-.00057 

295. 

-.00076 

-1196. 

-.00076 

-287. 

-.00065 

—  2  8  7  , 

-.00077 

356. 

.00047 

-283. 

-.00054 

283. 

-.00079 

-1217. 

-.00079 

-279. 

-.0006* 

-.0008? 

51  1  . 

.00056 

—  2 

-.00084 

21?. 

.00039 

-.00091 

-2856. 

-  ' 

-.00098 

-  "*  " 

•  >  1 . 

-.00198 

-27. 

- 

-10233, 

-.0019? 

10. 

,00209 

-6704. 

-.00106 

6, 

-668, 

-.00158 

-9190. 

-.00175 

-3e, 

.00068 

1626. 

.00188 

-6443, 

-.00169 

55. 

.0006? 

-706. 

-.00151 

-8494, 

-.00163 

-71. 

.00049 

1377. 

.00166 

-6181. 

-.00151 

103. 

,00048 

2012. 

-.00004 

-64  8. 

-.001*5 

-979. 

.00043 

-592. 

.00005 

-376. 

-.00121 

85  3. 

.  0 1/  0  3  8 

-2406, 

-.000)3 

-106. 

-.00097 

727, 

.00033 

463. 

.00005 

166. 

-.00073 

-601. 

CURVATup£S 

-.137E-0?  -. 

758E-C3  -, 

,  7  0  ?  F  -  0  2 

•  00038 

2602. 

.00369 

14036, 

.00359 

391, 

.00037 

-1864 . 

-.00327 

15503. 

.0035? 

-374. 

.  0  C  0  3  6 

23*9  . 

.00334 

13384. 

.00336 

340  . 

.00036 

—  1  P  ?9 . 

-.00309 

24347. 

.00328 

-323. 

.00034 

221*. 

.00317 

12940. 

.003?0 

305. 

.00032 

-1794. 

-.00291 

13190. 

.00305 

-271. 

.0003? 

19P1  . 

.00292 

12296. 

.00297 

264. 

.0003! 

-1770, 

-.00279 

12419. 

.00209 

-237. 

„ 0on?9 

17-'.  9. 

.00267 

11*43. 

.00274 

203. 

.10028 

1008. 

.0006? 

1945, 

.00266 

-1620, 

.00025 

-1632. 

-.00067 

1618, 

.00234 

14  46. 

.00021 

-1684. 

-.00060 

129?. 

•  0  0  2ft? 

1269. 

•  00ft]  p 

-2196. 

.00011 

966. 

.00169 

-1093. 

GM1? 


,ono6o 

,ono60 

.noose 

.00057 

.BOOST 

•nooss 

•  00054 
.00ft64 
,0M6? 
.00143 
,30164 
.00184 
.00204 


,  O  ft  o  6  0 
.00060 
.oonsft 
.none? 
.00057 
,00055 
.00054 


♦ 

>.00005 

.none? 
. ocooi 
.00007 
.onpi  o 
.00014 
,00020 
.00163 
.00142 
.001 5*1 
.00100 


.00075 
• .  0  fi  n  7  2 
,00065 
,0006? 
. 00069 
,0006? 
,  o  o  r.  a  9 

.0004* 

,000 39 
, C027G 
,00241 
.  C  C  ?  1  l 
.00152 
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Figure  (1).-  Composite  Cross  Beam  Finite  Element  Model, 
Drawn  Using  SPAR. 


Figure  (2).-  Undeformed  Composite  Cross  Beam  Finite  Element  Model 
With  Hidden  Lines  Removed,  Drawn  Using  HIDLIN. 


Deformed  Composite  Cross  Beam  Finite  Element  Model 
with  Hidden  Lines  Removed,  Drawn  Using  HIDLIN. 


Figure  (3) 


ite  Element  Model  of  Supersonic  Cruise  Aircraft 


Figure  (6) .  Deformed  Finite  Element  Model  of  Supersonic  Cruise  Aircraft,  Drawn 
Using  HIDLIN. 


drawn 

HIDDEN 

LINE 


Figure  (8) Composite  Cross  Beam  Finite  Element  Model  With  Most 
of  The  Hidden  Lines  Removed,  Drawn  Using  HIDLIN . 


CDC  SYSTEM  CONTROL  CARDS 


Figure  (9).-  Typical  Program  Setup. 


Figure  (10a).-  ax  Stress  Contours  of  Upper  Flange  of 
Composite  Cross  Beam  Structure. 


a 

Y 


Figure  (10b).-  cry  Stress  Contours  of  Upper  Flange  of 
Composite  Cross  Beam  Structure. 


Figure  (10c).-  Txy  Stress  Contours  of  Upper  Flange  of 
Composite  Cross  Beam  Structure. 


GEOMETRIC  NEUTRAL  AXIS 


Figure  (12) Program  SAT  Deck  Setup. 
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